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Abstract 

In this article, we study the no-boundary wave function in scalar-tensor gravity with various 
potentials for the non-minimally coupled scalar field. Our goal is to calculate probabilities for 
^ , the scalar field - and hence the effective gravitational coupling and cosmological constant - to 

m 

1/^ ■ take specific values. Most calculations are done in the minisuperspace approximation, and we 

^^ ■ use a saddle point approximation for the Euclidean action, which is then evaluated numerically. 

t~^ ■ We find that for potentials that have several minima, none of them is substantially preferred 

o: 

by the quantum mechanical probabilities. We argue that the same is true for the stable and 
the runaway solution in the case of a dilaton-type potential. Technically, this is due to the 
inclusion of quantum mechanical effects (fuzzy instantons). 

These results are in contrast to the often held view that vanishing gravitation or cosmological 
constants would be exponentially preferred in quantum cosmology, and they may be relevant 
to the cosmological constant problem and the dilaton stabilization problem. 
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1 Introduction 

One of the big challenges of contemporary physics is to obtain a working theory of quantum gravity. 
But even with such a theory in hand there could presumably remain questions of why the initial 
conditions of the gravitational field where what they were. Such questions may even pertain to con- 
stants of nature and the like, in a theory of everything, including gravity. It is therefore interesting 
to look for ways in which the initial conditions can be specified in a natural way. 

With this in mind, in the present article we study quantum Brans-Dicke theory 9,. In this 
theory, gravitational coupling and cosmological constant can be dynamical: With a non-trivial 
potential, one can typically find solutions of the classical equations in which the scalar $ carries out 
small oscillations around a minimum in the potential. Such a minimum can then be thought of as 
determining a pair of effective constants (G, A) . We note that a Brans-Dicke theory with potential 
can be obtained from string theory upon introduction of the dilaton field [151 HH] • The dilaton field 
is then related to coupling parameters of all interactions in string theory, and questions about the 
value of the gravitational coupling strength are thus connected to the dilaton stabilization problem^ 
i.e., the problem to explain why the field value of the dilaton of our universe is located in such a 
stable location |24) . This is a problem that has currently no commonly accepted answer. 

Upon quantization of the theory, each quantum state gives probabilities for the values of $, 
and these can, provided the state allows for a suitable classical interpretation, be translated into 
probabilities for effective constants G and A. Such quantum states are obviously not unique, but 
there are proposals to single out a unique natural state. Given such a proposal, one obtains an a 
priori probability distribution for G and A. 

For the canonical formulation of quantum gravity, the no-boundary (or Hartle-Hawking) wave 
function *hh [11 El [3] is such a proposal for a natural initial state of the universe!^ 

Probabilities for the values of coupling constants have been discussed in the context of quantum 
gravity for a long time. Early works (sec for example |25| ) rely on the fact that in the instanton 
approximation, the no-boundary wave-function is dominated by the Euclidean de-Sitter solution, 
whose action comes out proportional to 1/A, and argue that a vanishing cosmological constant is 
preferred by the resulting probabilities. A famous and detailed argument by Coleman [26^ has 

"5*0(1 oc-^^ (1) 



^It should be noted that there are obvious questions regarding how to interpret ampUtudes and probabilities 
calculated from a wave function of the universe such as ^hh- There is an extensive literature about and range of 
approaches to this question, but we will not address any of this here. Rather, we will take a pragmatic approach, by 
just assuming that the extreme cases (probabilities very large or small, similar or very different) do tell us something 
about the kind of space times that are excluded, or, respectively, described by a specific model. 



where A is to be interpreted as an effective constant that comprises contributions from a path inte- 
gral over wormhole solutions, and the exponential of (minus) the effective action is to be interpreted 
as a probability distribution for A. This result has also be interpreted as driving the gravitational 
coupling to zero [27] . 

Coleman's argument has been studied in the context of scalar-tensor theory [35] , with the same 
result for A, and the additional one that the whole theory would be driven toward the Einstein 
gravity limit. 

In the present work, we are more moderate and more ambitious at the same time. We do 
not consider the sum over topologies, and thus the coupling constants in our action will not receive 
corrections. On the other hand, as described above, the scalar determines the gravitational coupling, 
and introducing a potential, we obtain dynamical vacuum energy. Using the no-boundary wave 
function, we obtain probability distributions for both. We work in the minisuperspace setting 
[U m [3j |4] . This has been considered previously for Brans-Dicke theory, for example in [lOj [Til 
[T^ [T51 HM [T51 [TBI [T7] . In (TSl [17] various initial conditions for the wave function are studied in 
regards to the consequences for the effective cosmological constant, using a WKB approximation. 
The standard result of a minimal value for the cosmological constant is obtained for the tunneling 
proposal, but interestingly, not necessarily in the no-boundary proposal. What is new in our work 
is that we strive for a much better approximation to the no-boundary wave function. We follow 
[T] [8] in allowing fuzzy (i.e., complex- valued) instantons to contribute to the approximation, and 
we also use numerics to determine the solutions in the fully dynamical regime. The classicality 
condition of [7] [S] also plays a pivotal role. 

In Section[21 we introduce the general procedure to obtain an approximation to the no-boundary 
measure: Minisuperspace approximation, saddle point approximation of the no-boundary wave 
function, and imposition of classicality conditions. In Section [3] we apply the method to scalar- 
tensor gravity, and describe the numerical algorithm that we use to find the relevant stationary 
points of the Euclidean action. In Section [4] we report the results on the no-boundary measure 
for two types of potentials: quadratic potentials and double or multiple-well potential. Finally, we 
will apply our results to the dilaton stabilization problem. We have also included analytic results 
wherever we could. In Section [5] we summarize and discuss our findings. 



2 The no-boundary proposal in quantum cosmology 

2.1 Canonical quantum cosmology 

By reducing the full action to homogeneous and isotropic space-times, and upon choosing an ap- 
propriate time coordinate t the line element can be written 

ds^ = N^(t)dt^ - h.jdx'dx^ (2) 

where the positive definite metric h only depends on finitely many parameters q^{t), due to the 
symmetry assumptions. The corresponding action can typically be brought into the minisuperspace 
form 



S[N,q]^k / dXN 



A\ / ^ j„s 



1«^4^^ ail-"<^' 



(3) 



where q and N comprise the configuration variables, k is the coupling constant, q^ , I = 1, 2, ... n 
comprise the metric degrees of freedom, as well as the degrees of freedom of the matter fields, and 
1^ is a model-dependent effective potential. G is the DeWitt metric on minisuperspace. 

Let us describe the canonical formulation based on Equation ^. The canonical momenta 
conjugate to gare given hy pA = G Asq^ / N . N is non-dynamical. Its variation imposes an additional 
constraint 

C = \g^''paPb + V{q) - (4) 

on the variables {pA,? }■ The Hamiltonian H — NG vanishes identically on the constraint hy- 
persurface. While the kinematical phase space spanned by {pa, 9 } is 2n-dimensional, after going 
to the constraint hypersurface C = and identifying the gauge orbits generated by C, we end up 
with a 2n — 2 dimensional physical phase space. The system can be quantized a la Dirac, with 
commutation relations 

[q^,PB\=ihk8^l. (5) 

Upon choosing a suitable ordering, the constraint G can be turned into an operator, and physical 
states are required to satisfy the Wheeler-DeWitt equation G^ — 0. Solutions to this equation are 
typically not square integrable (zero is in the continuous spectrum of the constraint operator), but 
if G is selfadjoint, it determines a scalar product on physical states. In practice, it may be very 
difficult to obtain this scalar product explicitly. 

2.2 Classicality condition 

Let us assume for the discussion that follows, that the physical Hilbert space is a proper subspace 
of "H, and consider a physical state ^{q). Of particular interest are regions in configuration space 



over which ^ has a semi-classical form 

vl/(g)«A(g)e'^('), S{q)eR, (6) 

where the rate of change of S is much greater than that of A, 

|VM(g)| « |V/5((z)|, I = l,...n. (7) 

One way to see why a state of this form really describes almost classical behavior is the following: 
The Wigner function W[^] of ^ of a state satisfying Equations ([5]) and (O is approximately 

Wmq,p)^\A{q)\H{p-VS). (8) 

This shows that for q in this region, ^ determines a probability, and a momentum value p = VS* 
which has a high likelihood. Let us consider a surface S in minisuperspace that is intersected by 
each trajectory exactly once. Then the probability distribution on this surface in terms of q turns 
out to be (cf. ex. [7]) 

p{q)^\A{q)\^ n-VS = n-J (9) 

where n is the normal to S in minisuperspace and J is the conserved "Klein-Gordon" current 

ih ^ 
J = **V*. (10) 

2.3 No-Boundary wave function and steepest descent approximation 

In the minisuperspace approximation, the no-boundary wave function takes the form 

*HH(g) = / pg(•)e-^-l«(■'l/^ (11) 

where NB((7) are symmetric Euclidean space-time histories q{-) with a single boundary, the geometry 
of which is described by q. We assume that the parametrization of the histories, or equivalently, 
the lapse function, has been fixed appropriately. But even after gauge-fixing, the integral is not 
convergent because the Euclidean action is in general not positive. One way to possibly cure this 
divergence is to regard the path integral in Equation (fTTj) as a contour integral, and deform the 
path away from real superspace, into the space of complex symmetric metrics. Despite the huge 
simplification due to the minisuperspace approximation, the path integral (Equation (TT|)) is hard 
to treat exactly. Therefore one uses the steepest descent approximation. One obtainf|^ 

*HH(g) « Y. ^(9cxt)e-^-[«=-l/'' (12) 



^Note that, here and in the following, Se [q] will sometimes be considered a function of the endpoints of the history 
(?(•). We hope that this does not lead to confusion. 



where Qcxt is an extremizing Euclidean history satisfying the appropriate (no-)boundary conditions, 
and P is given by 1 plus higher order terms which depend on functional derivatives of the action. 
We will neglect the higher orders of P in what follows. For a single extremum (in the application 
below we will see that the boundary conditions allow at most two extrema) we then havqf 

Comparing with Equation ^ , we see that ^hh is of scmiclassical form at q if 

IV/^HcSeI « |V/3m^E|. (14) 

Let us finish by discussing a subtle point that is nevertheless important in applications: we have 
said earlier that to give meaning to the path integral, we have to gauge-fix the lapse, or equivalently, 
the time parameter used in the action. A priori, the time parameter is real, but it is very convenient 
to also consider complex time parameters. In fact, if the Lagrange function and the history are 
analytic in suitable regions, the action can be viewed as a contour integral in the complex time 
plane, and hence depends on the integration contour only through its endpoints. Following [6l[7] we 
will make use of the freedom to chose the contour, below. Since we have written the no-boundary 
wave function in terms of the Euclidean action, we will refer to a section of a contour that is parallel 
to the real axis as Euclidean, to one that is parallel to the imaginary axis as Lorentzian. 

3 The no-boundary measure in scalar-tensor gravity: For- 
malism 

In this section, we investigate the no-boundary measure for scalar-tensor gravity. Scalar-tensor 
gravity is of interest not only in itself, but also as a dilaton gravity limit of string theory. For us, 
it is interesting because it has dynamic gravitational coupling and vacuum energy. Quantization of 
scalar-tensor gravity has certainly been considered before. In the context of quantum cosmology 
see for example [ini HH [HI HSl HH HSl [161 [12 and our discussion in the introduction. For non- 
perturbative quantization, see for example [501 [SI [121 [23] ■ 

The first and simplest scalar tensor theory is that of Brans and Dicke [S] , 

Here, cj is a dimensionless coupling parameter, and Einstein gravity is restored in the w — >■ oo 
limit. We will refer to more general actions involving a non-minimally coupled scalar as scalar- 
tensor theories, whereas we reserve the name Brans-Dicke theory to the theory with the above 



'In this paper, IHe denotes the real part and 3m denotes the imaginary part. 



action, possibly with the inclusion of a potential V{^). Matching current observations against the 
Brans-Dicke theory without a potential requires w greater than ^ 40,000 [29 . Small values for w 
are however admissible for non-trivial potentials |30j . and can be found in string-inspired models 

[Mils]- 

We briefly comment on possible scenarios that involve Brans-Dicke theory. The most important 
one is dilaton gravity obtained from string theory. The effective action of string theory has the 
form [S]: 

S = ^ J d'^+'xV^e-f {R + (V^)2) , (16) 

with d the space dimension. As the length scale of string units, R the Ricci scalar, and cf) the 
dilaton field. The field redefinition $/(87rGd+i) = exp(— 0)/Af^^ turns the above action into that 
of Brans-Dicke theory with oj = —1. 

In the first model of Randall and Sundrum [32 , two branes have been employed to account for 
the hierarchy problem. Due to the warp factor between two branes, one obtains a positive tension 
brane and a negative tension brane in the anti de Sitter space background. According to Garriga 
and Tanaka [33j each brane can be described by Brans-Dicke theory in the weak field limit with 

3 



OJ = 



I (e±^/' - l) , (17) 



where s is the location of the negative tension brane along the fifth dimension, I — ^— 6/A is the 
length scale of the anti de Sitter space, and the sign ± denotes the sign of the tension. To explain 
the hierarchy problem, we require s/l ^ 35. We then obtain a sufficiently large value of w on the 
positive tension brane while uj > —3/2 on the negative tension brane |331I18| . In principle, however, 
s/l can be chosen arbitrarily, and hence one may infer that various w near —3/2 may be allowed by 
models of the brane world scenarios. 

Finally, it is well known that /(i?)-gravity (for a recent review see [M]) can be reformulated in 
terms of a Brans Dickc field. The action is 

1 



S — dx \J—g 



, „ f{R)+ -Cmatter 
iDTT 

Introducing an auxiliary field ip, we change the gravity sector to 

1 



(18) 



^gravity = J^ J dx^ V~9 [f {^) + f'WiR - ^ (W) 

with constraint t/j = R. If we define a new field $ by <& = f'{tp), we obtain 

^gravity ^^J dx^ V^ [*^ - V^ (20) 

where V{^) — ~-f{ip) + tpf'{ip). This is exactly the oj = limit of the Brans-Dicke theory. 



For our purposes, we begin from the action (Equation ([T5|)) by adding a potential V{^). In this 
article, we focus on the two types of potentials: quadratic potentials 

y($) = iAf2($ _ 1)2 = 1a/2s2, (21) 

where S = <I> — 1, and double- well potentials, 

V{^)^^'(£^d^ + Vo\ (22) 

where we use an effective force function F(<i>) given bj|f 

F($) EE $y' ($) - 2y($) = yl ($ - $a) ($ - $b) f $ - ( ^2^L^ + ^ ) ) • (24) 



A is a positive constant, $a and <i>b denote the field values of each vacuum, and 5 is a free parameter 
that determines the location of the bump of the potential. For convenience, we choose <i>a = 1 and 
let ^(<&a) — Vq be true vacuum value. It is also useful to introduce the effective potential 

F(l')d$=/ ($y' ($) - 2y(l>)) d$, (25) 

which is the one showing up in the field equation, V^^ = U' /{3 + 2cj). 

In this article, we will not explicitly perform calculations for a dilaton type potential resulting 
from the field redefinition in Equation (J16p . but conduct a qualitative discussion of the stability 
issue based on an approximation by a double well potential. For a quantitative study of the dilaton 
type potentials see [55] . 

We write the symmetric line elements aqj 

dslj^ = TN^it)dt^ + p^{t){dx^ + sin^ x{dO^ + sin^ Od^p^)), (26) 

and obtain the symmetry reduced actions 

S'l,e[/0, $, iV] = I I dtN {±Qp^^pN-'^ ± 6$pp2Ar-2 

-6$p T ujp^^^^^'^N-^ + pV($) 



(27) 



To obtain the above result, we have added a suitable total derivative to keep the action differentiable 
in the presence of boundaries, dt should be changed to d-q for Euclidean case. 



*Note that the potential in the Einstein frame V is I19| 



V(^) = ^ 



^3 d* + Vo. (23) 

The dynamical field <E> in that frame is related to <J> by <& = exp 'i>[(167r)/(2Lj + 3)]"'^'^. Therefore, the only effect is to 

stretch the potential along the field direction, and this does not affect the vacuum energy of each field value. 
^For some more detail regarding the following calculations see for example |36| . 



For the approximation of the no-boundary wave function, the on-sheh Euchdean action and 
both, Euclidean and Lorentzian equations of motion are needed. Using the Euclidean Einstein 
equation for p and field equation for </> (with A^ = 1), one finds 

Se^\J dr^{p^V-Qp^). (28) 

The equations of motion (with iV = 1) for the Euclidean time rj and the Lorentzian time t are 
as follows: 

$ = -3^$ ± -i— ($V^' - 2y) , (29) 

p ZUJ + 6 

±1-P' <&' ^ pV _ 2up V 

where the upper signs are for Euclidean and the lower signs are for Lorentzian. Note that all 
functions in the equations are complex in general [p — p^'^ +ip-^'" and $ = ^^'"^ +i$-''"). Therefore, 
effectively, there are four functions (p^% p^"', $^% $•'"') and we need eight initial conditions to fix 
a solution. 

As we have already sketched the general procedure in Section 12.11 we can be brief here. The 
symmetry reduced Lorentzian action, Equation (j27p . is of the form Equation ([3]), where q = (p, <&) 
and N comprise the configuration variables, the coupling constant takes the value k — tt/8, the 
effective potential is given by V^ = 6p<l> — p^V, and 

Gpp — 12p<f>, G^p = Gp<^ = 6p , G$$ — —2a;— (31) 

are the components of the DeWitt metric on minisuperspace. Thus the system is equivalent to a 
relativistic particle propagating in a two-dimensional space with metric Gab under the influence of 
a potential V. The determinant of the metric is — 12(2a; + 3)p'*, so Gab is Lorentzian for uj > —3/2 
and Euclidean for a; < —3/2. 

The detailed expressions for the momenta and the Hamilton constraint C do not concern us 
here, but let us make some remarks on the reduced phase space. As with all reparametrization 
invariant systems, time evolution in gauge evolution. While the kinematical phase space spanned 
by {pA,q } is 4-dimensional, after going to the constraint hypersurface C — and identifying the 
gauge orbits generated by C, we end up with a two-dimensional physical phase space. Points in 
this space are simply the trajectories $(jo) that can be obtained by solving Equation (|3]) to obtain 
q'^{X) and then eliminate A which is possible at least locally. 

Let us also give a parametrization of the space of trajectories ^{p). Let us assume that the 
particle is not a tachyon - the parametrization can also be easily adapted to the more general case, 
but details then depend on the potential. We fix an initial value slice p ^ pa = const. Note that 



this slice is spatial with respect to G. Then each solution of the dynamics will intersect this slice 
once, and we obtain $(/Oo) find the two momenta p/i(po) £^t the intersection point. We note that 
reparametrization of the solution changes the momenta by a constant factor, so we can chose as 
parameters $(/9o) and the ratio Pi{pq)/p2{po)- 

3.1 Steepest descent approximation 

We remind the reader that to obtain the steepest descent approximation of the no-boundary wave 
function (see Section [2. 3|) at a superspace point p, $, we need to determine complex solutions of the 
equations of motion that have p and $ as boundary values and otherwise fulfill the no-boundary 
condition. To specify the boundary value problem fully, and to compute the action of such a 
solution, we need to fix a contour in the complex time plane. While this choice is in principle largely 
arbitrary, there are practical reasons to chose one contour over the other. To formulate the no- 
boundary condition, it is necessary to start the contour as Euclidean. Moreover, the interpretation 
of the results is done most easily if the end of the contour is Lorentzian. 

The contour we chose for the practical evaluation consists of two components: An Euclidean 
component from rj ^ to rj = X and a Lorentzian component from t] = X to j] = X + iY (or from 
t = to t — Y with rj = X + it). To implement the no-boundary condition, we require that 

p(Of ^ - p(0)3- = 0, P(0)«^ = 1, piO?"' = 0, (32) 

4(0)^^ = $(0)^"^ = 0. (33) 

At the turning point j] = X, we have to match the Euclidean part of solutions (so called, p, p, 
$, $) to the Lorentzian part of solutions (so called, p, p, $, $) of Equations ([^ and ([5(11 via the 
matching conditions: 

p{t = Q) = p{rj = X), p{t^O)=zp{f^^X), (34) 

$(i = 0) = $(r/ = X), $(t = 0) =i$(?7 = X), (35) 

since d/dt = id/ drj. 

With a view to the calculation of the no-boundary wave function, $ (/?,$) we note that we 
have 6 real conditions at 77 = 0, which leaves two real initial values to be specified. Additionally, 
our contour has two parameters {X,Y). This makes 4 real parameters that we can shoot for at 
the boundary, i.e., precisely the number needed to tune the arguments of the no-boundary wave 
function. We are, however, only interested in regions of configuration space on which the wave 
function has semiclassical form. This is what we discuss next. 



10 



3.2 Classicality condition 

We are not interested in the value of the no-boundary wave function for arbitrary arguments, but 
only in regions of configuration space on which the wave function has semiclassical form Equa- 
tions © and ([7]). We note that the latter condition actually has two components, of which one can 
be checked easily in our setup. The easy component is the one in the direction of time evolution: Let 
q{t) = (a(i), $(i)) one of the no-boundary histories, evaluated on part of the Lorentzian contour. 
With dq — dq/dt\tg, we find 



Sq^WASB{q{to)) = ^ 



C{q{t'),q{t'))dt' = £{q{to),q{to)), (36) 



where C is the Lagrange function. Thus, given a history q(t), the component of Va5'e in the 
direction of 5q is easy to determine. We also note that, by definition, in a classical region q and 
6q have to be real. In what follows, and in slight abuse of terminology, we will call a no-boundary 
solution of the equations of motion which satisfies the tangential component of classicality, with Sq 
real, a classical history. 

The orthogonal component of the classicality condition is harder to check, as it involves finding 
no-boundary solutions of the equations of motion with nearby boundary conditions. 

For the numerical evaluation, we split the initial value of $ as 

$(0) = $oe*^ (37) 

with $0 and 9 real. We fix the latter, and then numerically scan the {0, X) space for solutions that 

1. become classical at late times, and 

2. have p''"', p , $ , and $ approach zero in the large t limit. 

The parameter Y is held fixed (and sufficiently large), as it simply corresponds to time evolution 
in Lorentzian time. 

3.3 Probabilities 

We will denote the space of histories satisfying the no-boundary condition by i^O We eventually 
want to calculate the probability of histories satisfying certain (boundary-) conditions. Let us say 
that we are interested in histories satisfying a certain condition A. Then we can define the subset 



®Note that histories which differ by time reparametrization are considered the same. Notice also that while these 
histories are regular in Euclidean time as per the no-boundary conditions, they may well have singularities along the 
real time axis. 



11 



by 

^A = {h E S) \ /i has property A} (38) 

of Sj. Given a set of classical histories ^a, in principle we have 

Pa= f \'^Mh,x)?n■VSV^l{h,x). (39) 

JQa 

The integration is over a subset Qa of a spatial slice with normal n in superspace. Qa is the 
set of points on this slice such that the wave function satisfies the classicality condition in a way 
compatible with the condition A. S was defined in Equation ([B]). fi is a certain measure which can 
be obtained in principle from the inner product on the space of solutions to the Wheeler-DeWitt 
equation and the slice. But it is very difficult to obtain in practice. Using minisuperspace and 
steepest descent approximation, using $0 as parameter on the slice, and ignoring details of the 
measure as well as the variation of n ■ VS*, by assuming both are constant over the space of histories 

Pa~^ f |exp(-5ENJ)prf$o = I / e-2«c5K[/..J^$^ (40) 

'^ JQa ^ JQa 

where Z is some normalization constant and hi^^ is a history that initially has scalar field modulus 
equal to $o- 

If we have two conditions A, B, where A is an initial and B a final condition, we will use the 
notation ^a^b for ^a H S)b, and Pa^b for the corresponding probability. 

If S'e[^$o] i^ slowly varying, and we compare two probabilities, then from Equation (|40p we get 

-^ « exp (-25e[%J + 25e[%J) ^^^ ,. , (41) 

where <i>i is a suitable initial condition such that the history fulfills Ax, and $2 the same for ^2- 

3.4 Searching algorithm 

We have to find initial conditions, the initial phase angle 6, and the turning point X, for a given 
initial field amplitude $0 to satisfy the classicality condition. To find classical histories for a given 
$0j we formulate an optimization problem and solve it using the idea of the generic algorithm. 
To realize this algorithm, first we define the objective Junction that quantifies classicality. Second, 
using the optimization algorithm, wc list and choose the best candidates of initial conditions using 
the scores of the objective function. Third, we inversely check whether the solutions of the initial 
conditions really have classical properties what we required. 

We first define the standard objective function F^^. To this end, note that the tangential 
component of classicality can be reformulated as 

V(Va5e)^^ _ C^'iqA) 



Sq^VaSe 
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« 1, (42) 



as long as Sq^ is real (see the discussion around Equation ([5S)) ). Evaluating the left-hand side of 
Equation (1421) in a single point given by a numerically determined function is not a very stable 
procedure. Thus, in practice we average it over some time interval, to define the objective function. 






dt. (43) 



Ti and T2 define the Lorentzian time interval where the classicality will be tested, and need to be 
chosen sufficiently large. The reality of 5q^ is checked separately for minima of F. Altogether we 
obtain a criterion that is, for asymptotically large T, strictly equivalent to one component of clas- 
sicality. Other definitions for objective functions that single out classical histories are undoubtedly 
possible. As said sketched above, ours was chosen for numerical conveniencql ■ 

Then, we can define the relevant optimization problem: for given $0 a-nd other constraints, which 
values of initial parameters {0,X) minimizes F^^l Of course, we can find the optimal solution by 
searching all possible values of {9,X). However, it takes too much time and computation power, 
and hence we need a better optimization algorithm. 

In this paper, we use a simple version of the generic algorithm. 

Initialization: We first generate a set of initial condition pairs of a number N: {{Oi^Xi)} {i — 
1, ...,N), where we choose 9 and X randomly. This set of {Oi,Xi) constitutes the first gener- 
ation of candidate solutions. 

Elite: Then, we calculate F^^ for each (0^, Xi) by solving equations of motions numerically. Among 
the parameters {9i,Xi), the one that has smallest F,j(, can be regarded as the most classical 
solution in the first generation. As a subset of candidates, we select the Ug < N number of 
initial conditions that have smallest F^g . We call them elites of the first generation. 

Cross-over: We cross-over the elite parameters to generate the Uc [n^ + nc < N) number of new 
parameters. That is, we pick two parameters (0^, Xi) and {6j, Xj) for arbitrary chosen i and j 
from the elites and define a new parameter {{9i + 9j)/2, [Xi + Xj)/2); we repeat this process 
Uc times. They are adopted to examine the parameter space around the elite group more 
detail. 

Mutation: Finally, we generate the Um {n^ + nc + n.m ~ N) number of totally new parameters, 
so-called mutations, by choosing 9 and X randomly. They ensure us that potential candidates 
of parameters are not lost systematically. 



'^We have run checks with other possible definitions for the objective functions. We found that the results generally 
agree, although some differences are observed in the performance of the numerical search algorithm, in particular in 
some more extreme regimes. 
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Evolution: The set of elites, cross-overs, and mutations constitutes the second generation. Using 
the second generation of initial conditions, we calculate the objective function F^^ . Then we 
can choose new elites, cross-overs among the new elites, new mutations, and define the third 
generation. We repeat this process to evolve generations until the values 9 and X of the elite 
group reaches to a steady state. 

This generic algorithm finds the optimal parameter as the number of generations becomes sufficiently 
large. However, it is possible that there is no classical solution for a given $o- Therefore, we have 
to be careful whether the optimized value is really the classical solution or not. 

There are two main drawbacks of the generic algorithm in searching classical solutions. One is 
that it can converge to a local minimum of F^^ rather than the global minimum, or there can be 
two or more classical histories for a given $o- To resolve this problem, we divided the searching 
region of (9, X) into several pieces and found the optimal parameters separately. In this way, one 
can reduce the possibility of being captured by a local minimum. The other drawback is that the 
result can sensitively depend on simulation parameters: N, n^, ric, rim, and the definition of the 
objective function. N, Ug, ric and rim should be sufficiently large until the result is not sensitively 
depend on the choice of them. In our choice of -F$q , we have to choose a proper time interval [Ti , T2\ 
where we will regard that a history becomes classicalized around the time. Therefore, we have to 
carefully choose the time interval case by case and check whether the results arc not sensitively 
depend on the choice of parameters. 

4 The no-boundary measure in scalar-tensor gravity: Re- 
sults 

In this section, we study the no-boundary measure for two types of potentials: quadratic potentials 
and the double-wcll/multiplc-well potentials. The former is useful to obtain clues abut the physics 
the near a generic local minimum, the latter to compare the probabilities for various values of 
gravitational couplings In the case of the quadratic potential, we also have some analytic results 
that we can compare to the numerical ones. Our results on these types of potentials will give some 
intuition for the dilaton stabilization problem. First, we will make some analytic considerations, 
then we will turn to the numerical investigation. 
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4.1 Analytic considerations regarding classical histories 

For a first estimate, we regard that the potential is approximated near $ = 1 by 

F($) ~ -Af2($- 1)2 + C'(($-1)^). (44) 

After we re-define the field S = $ — 1 , the on-shell Euclidean action becomes 

SE^lJdij(^lp'M^E^-6pE~6py (45) 

Therefore, compared to the Einstein case, we have an effective new term proportional to — pS. The 
field equation is 

E = -3-E + m^E, (46) 

P 

where we have defined m^ = M^/(2w + 3). Thus, for a quadratic potential the Brans-Dicke field 
equation is the same as that for a standard scalar coupled to Einstein gravity. Thus, as it was 
approximated by |5j, we can compare the Euclidean action to the slowly rolling and almost static 
field limit. 

In the static case, where $ = $o = const, with $o determined by $oV^'($o) ~ 2y((f>o) = 0, we 
have the exact solution, 

$(?7) = $o, p{7]) ^jsin Arj, (47) 



where 



A^ 



1 



V = ^V. (48) 



_4a; + 6 6a; + 9 

To satisfy no-boundary and classicality conditions, the turning point of the contour should be at 
X — At:/2. Then the real part of the action does not grow after the turning point and the action 
becomes 



X Q^<fi2 



5e = J y dTj (p3y($o) - 6p<I>o) = - 



37r$, 







(49) 



n*o)' 

Therefore, for small S limit, the Euclidean action negatively increases and hence $ '^ 1 will be 
preferred. 

Let us go even a little further: We have found one exact solution fulfilling classicality, and it is 
a de Sitter-like solution, i.e. trigonometric/exponential dependence of the scale factor in Euclidean 
respectively Lorentzian time. Thus it is not unreasonable to guess that other solutions contributing 
to the classical regime are perturbed de Sitter like solutions. To obtain those, we need a field 
configuration allowing an inflating space-time. The well-known mechanisms are slow-roll inflation 
and false vacuum inflation. We are therefore led to hypothesize that histories contributing to the 
classical regime are allowed under the following conditions: 
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1. A part of the potential satisfies the slow- roll conditions and the history experiences the region 
of the potential (thus mimicking slow- roll inflation). 

2. The field becomes slowed-down and approaches a local maximum or minimum to experience 
an effectively de Sitter space (thus mimicking false vacuum inflation). 

The slow-roll conditions for Einstein gravity are 

2 




(50) 



In terms of the Brans-Dicke potential, 



167r \ /„ ^V'V . f 16tt \ ( , „^V' ^^V" 



From this consideration, one is led to the following hypotheses: First, assuming the local extrema 
to be obtained for l$| ^ C(l), both of the slow-roll conditions are difficult to satisfy, unless there is 
an inflection pointO Second, assuming V ^ $" for large $, slow-roll conditions hold only if n = 2. 
Therefore, if the potential shape is "runaway" for large $, i.e., n < 0, classicalization via slow-roll 
inflation may not happen. Certainly, in both these cases, classical histories may still appear via 
false vacuum inflation. 

For the quadratic potential, the conditions are 

Therefore, both of the conditions hold if and only if in the $ ^ I limit. This suggests that there are 
no classical histories that start near the bottom of the potential $o ^ 1- Such behavior has been 
observed numerically in the Einstein case [7] . It also suggests that the range of initial conditions that 
allow for classical histories gets larger as uj increases. Both of these conclusions from our analytical 
considerations are borne out very well in the numerical results that we will present below. 

4.2 Numerical results on the quadratic potential 

We have numerically determined the classical histories contributing to the no-boundary condition in 
the case of the quadratic potential. Our algorithm converges well and finds a region of initial values 
for which there are classical histories. In this region one solution per initial value $o — \^{v = 0)| is 
found. This is in contrast to the Einstein case, where there are two solutions due to the symmetry 
of the potential, under i— >■ — (/>. Figure [1] shows an example of a classical solution in Euclidean and 
Lorentzian time, respectively. In Euclidean time, the real part of the scale factor p increases as a 



*Maybe an exceptional case can happen if the potential is ~ <I'" + A and A is larger than 1; however, a cosmological 
constant larger than 1, is unrealistic for our universe and we are not interested such a case. 
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Figure 1: The classical solution p^', p^^, $^% and $•''" as functions of rj and i, for w = — 1, 
M = 0.75, and $o = H- In this case, the search algorithm determined the turning time to be 
X = 2.0734 and the initial phase angle to be 6* = 5.68017. 
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sine function, which is typical of a de Sitter-like space in the Euclidean signature. The real part of 
the field $ slowly increases, since for Euclidean signature, the potential is effectively inverted. After 
the turning point, the real part of the scale factor increases exponentially as is typical when there 
is positive vacuum energy. $ rolls down to the equilibrium near $ = 1, as expected. Other classical 
histories show similar characteristics, as we had already anticipated by analytical arguments. 

We also note that the solution becomes real in the large t limit, as it should; It is instructive 
to compare the real part and the imaginary part (top and middle of Figure [2|). After the turning 
point (marked by the red circles) , the imaginary parts quickly decrease to zero and are in particular 
extremely small as compared to the real parts. 

Also the behavior of the action (bottom of Figure [2]) is consistent with our expectations. During 
the Euclidean time, the real part of the Euclidean action increases (to the negative direction) and 
variation of the imaginary part is negligible. However, after the turning point, the real part of 
the Euclidean action is almost constant, while the imaginary part of the action increases (to the 
negative direction). This shows that for large t, the variation of the imaginary part of the action is 
much bigger than that of the real part, and hence the classicality condition holds for the solution. 
The value of the Euclidean action at late time t can then be used to determine the probability of 
such a history: P ^ exp(— 2S'e). 

In Figure [21 we plot the Euclidean action as a function of |S(0)| for several values of cu. The 
behavior is as anticipated in the previous section: as $0 increases, probability decreases. For a 
given uj, there is a critical $0 so that there are no classical histories with initial value less than 
the critical one, which confirms our discussion using the slow roll conditions. As we decrease uj, 
classical histories can be observed in the small $0 region. This is related to the effective mass 
m^ = M^/(2aj + 3). If the effective mass increases, or cj > —1.5, then the field rolls more quickly 
than in the small effective mass cases, and hence it will be difficult to see a classical history near 
the minimum. 

4.3 Double-well potential 

The double- well potential has two minima for $; this is very important for the present work since 
it gives a simple model that allows two different universes with different gravitational couplings. 
After we gain intuition on the no-boundary measure in this case, we can apply it to the dilaton 
potential, by shifting one minimum to infinity. 
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Figure 2: Plots for p^'-p^"", $«^-$3"i, and S^'-Sl"" for w = -1, M ^ 0.75, and $o = H- Red 
circles are the turning point. 
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Figure 3: Euclidean action Se for quadratic potentials with w = 0, —1, —1.4. 
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Slow-roll Conditions 





Figure 4: A typical double- well potential in the Jordan frame (upper left), Einstein frame (upper 
right), and slow- roll conditions (Equation ([FT]) , the second of them positive definite) of the potential 
(lower left), for A = 1, $a = 1, ^b = 2, 5 = -0.05, and Vq = 0.0001. Lower right: The regions near 
the local minima (black regions) do not attract trajectories during the Euclidean evolution. The 
local maximum (green region) is an attractor during Euclidean time. 
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4.3.1 Analytic considerations 

Let us consider a simple case 

F($) = $y'($)-2V"($) (53) 

= Ai<i>^<P,)i^-<^^)U-f^^±^ + 6\Y (54) 

and ^ = 1, $a = 1, $b = 2, (5 = -0.05, and Vq = 0.0001. Wc plot the potential in the Jordan 
frame V{^) (upper left of Figure S]) , potential in the Einstein frame V{^) (upper right of Figure S]) , 
and both of slow- roll conditions (lower left of Figure S]). One can see that the slow-roll conditions 
never become sufficiently small at the same time. Perhaps, if a saddle point and an inflection point 
coincide in the Einstein frame, it may be possible to obtain a region where both slow-roll conditions 
hold. However, in this case the solution would not be stable and hence we will not consider this 
possibility further. 

However, even though slow-roll conditions do not hold, if initial conditions of the field are finely 
tuned so that if the field is slowed down and approaches a local extremum, there may be a possibility 
to see a classical history. One trivial example of this is the solution in Equation (H71) . Therefore, 
one may guess that such a real solution can happen near the solution of ^V'{^) — 2F($) — (local 
minimum or local maximum of the potential in the Einstein frame) . 

However, if the solution is a minimum, this can not work. If the field does not start exactly at 
the minimum, it will roll to the upper region of the potential during the Euclidean time evolution. 
Only a local maximum is an attractor during Euclidean time (lower right of Figure |4]), and may 
thus provide for a continuum of classical histories via false vacuum inflation. Thus there is some 
hope to see classical solutions 

1. precisely at a local minimum, 

2. precisely at a local maximum, and 

3. for initial conditions near the local maximum. 

Now let us define subsets of all histories f) by imposing initial conditions: S)ai, ^A2i ^-nd ^Aai 
such that 

JOAi = {/l|$0 = $m}, (55) 

^A, = {h\^o = *m}, (56) 

^A, = {/i| $o>*M-A$(i) or $o<*M + A$(2)}^ (57) 

where $m is a local minimum, $m is a local maximum, and A^^^^ A$(^) > are bounds to 
allow classical histories around the local maximum. We know that around the local minimum or 
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maximum, the Euclidean action will be approximately Equation (j49p , and hence 

Pa.. " m(^A3)' ^^ 

where /i is given by the Lebesgue in the space of initial conditions. Therefore probabilities will be 
dominated by Pa.i since sets JOi and ^2 are measure zero sets. 

Now let us focus on the Case {in). Then there may be two possibilities: after the history has 
become real- valued, it can roll down to left or right in the double-well potential. Let us impose the 
final conditions Bl and Bij and define subsets S^Bi^ and S^Br by 

^B^ = {h I $(A = 1) = left side}, (60) 

S^B^ = {h I $(A == 1) = right side}, (61) 

where A = 1 means a sufficiently large time along the Lorentzian direction. If the potential near 
the local maximum is symmetric, then for given $0i both of histories will be allowed. In other 
words, after the field is slowed-down at the top of the hill, if the velocity of the field is almost zero, 
then there will be no principle to push the field left or right. Therefore, as long as the potential is 
approximately symmetric near the local maximum, if a left-rolling solution is allowed, then there 
will be a right-rolling solution, too. In other words, 

and there is no exponential contribution to determine left or right. 

4.3.2 Numerical confirmations 

We have numerically confirmed these theoretical assertions. For convenience, we choose parameters 
to almost symmetric near the local maximum: v4 = 1, $a = 1, 'I'b = 2, (5 = —0.05, and Vb = 0.0001. 
This potential has local minima around $ = 1 and $ = 2 and the local maximum is <& ~ 1.45. 
Figure [5] is an example of the left-rolling solutions and Figure [5] is an example of the right-rolling 
solutions. In both cases, we fixed $0 — 1.4. For the right- rolling case, since the final state has 
sufficient vacuum energy, one can see the exponentially increasing p; therefore, future evolutions 
of left- rolling and right- rolling cases are quite different. This implies that for a given initial field 
amplitude $07 there are two physically different solutions (with different and X). Therefore, after 
we fix the both of initial and final conditions, $0 points out a unique history. In both cases, along 
the Euclidean time, the imaginary parts are definitely suppressed. 



22 





-real part 

- imaginary part 



/ 






reaJpart 


\ 







1.50- 
1.25. 

1.00- 
0.75. 
0.50- 
0.25. 
0.00- 




— real part 

— imaginary part 



Figure 5: An example of left-rolling solution (0^% p-'''", $^', and $-'™ as functions of rj and t, for 
A=l, $a = 1, $b = 2, (5 = -0.05, Vo = 0.0001, and $o = 1-4. 
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Figure 6: An example of right-rolling solution p^*^, p-'"', $^% and $*'™ as functions of r\ and t, for 
A = 1, $a = 1, *b = 2, (5 = -0.05, Fo = 0.0001, and $o = 1-4. 
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Figure 7: Euclidean action ^e for double-well potential: A = 1, $a = 1, $b = 2, (5 = —0.05, and 
Vo = 0.0001. 

As we vary ^q, we can classify left- rolling solutions and right- rolling solutions and estimate 
Euclidean actions. The allowed region for left-rolling solutions is $o ^ 1-8 and the allowed region 
for right-rolling solutions is approximately $o ^ 1-4 (Figure [T]). Moreover, there is no meaningful 
difference on actions between various $o and left or right-rolling processes. The values are approxi- 
mately — 1349, where — 37rl.45^/V^(1.45) ~ —1348.9. This confirms our assertions on the probability 
(up to the overall normalization factor Z): 



Pa3^b ^ exp 



and 



Pa, 






rf$o 



o(l), 



(63) 



(64) 



PA,-^Bn IhenA.^B^ '^^" 

where $m is a local maximum in the Einstein frame. 

In summary, we find in all the potentials that we have investigated numerically, that if there 
is a local maximum between two local minima, then the probability is determined by the local 
maximum. 
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4.3.3 Generalization to the multiple-well potential 

Now let us generalize for multiple-well cases. If we see a triple-well potential, then there are 
two places with non-zero measure where classical solutions can be obtained: let us call ^Ai and 
^A2 (Figure IS]). Then, there can be three final conditions: ends at the first minimum .^i, ends 
at the second minimum f)2, and ends at the third minimum f)^. Histories for Sji should begin 
from Ai, and each history contributes the probability approximately exp (67r$^ /V{^Ai))- Histo- 
ries for Sj2 can begin from Ai or A2 and each history contributes the probability approximately 
exp {Qtt^\ /V{^Ai)) or exp {6tt^\ /V{^a2))- Finally, histories for ^3 can begin from A2 and each 
histories contribute approximately exp {6n^\ /V{^a2))- Now we write probabilities by (we ignore 
the prefactor l/Z), 

P^, ~ exp— -^ / d$o, (65) 



P,,2 ^ exp— -^/ d^o + e^p-—^ d$o, (66) 

V^(*Ai) AeflAi^Bi^ V(^A2) Jh&S^A2^B2i^ 

P^3 ^ exp— -^ / dc&o. (67) 

y^^A^) Jhe^A2^B2n 

This argument can be generalized for other multiple-well potentials. 

4.4 Stabilization problem: dilaton-type potential 

Let us think a typical dilaton-type potential (Figure |9l). In many models of dilaton potential, 
there is an unstable direction in the large $ limit. Let us try to consider this as an extreme limit 
of a double- well potential, as in Figure IHl Let us call the minima $mi and ^m^ . As we have 
seen previously, histories in which $ is precisely located in one of the minima will lead to classical 
points of the no-boundary wave function with modulus squared proportional to exp67r$^/T^(<I>m). 
Therefore, if we compare just these two points, the right minimum '^m^ is much more preferred. 
We can think that a dilaton potential is an extreme limit of '^m.^ — > cx). In this point of view, 
Euclidean quantum cosmology seems to say that the dilaton field should be destabilized and all 
coupling constants of nature should be zero. 

However, as we discussed in the previous section, other histories also contribute classical points, 
in particular the fuzzy instantons. For the potential in Figure [9l it is also difficult to find a region 
that allows slow-roll conditions to be satisfied, except infiection points. Thus there can be classical 
solutions around the local minima or maxima; those starting near but not precisely at the local 
minimum will evolve towards the maximum and obtain the main contribution to the Euclidean 
action there. Furthermore there will be continuously many classical histories starting around the 
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Figure 9: A typical dilaton-type potential and an extension for a thought experiment. 
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local maximuni of the potential. Therefore, as in Equations (j58l) and (j59p . the local minima will be 
excluded. 

For obvious practical reasons, we cannot directly work with the $ = cx) solution itself. Numer- 
ically, we only see runaway solutions. However, as we numerically observed, all runaway solutions 
should experience false vacuum inflation around the local maximum, and hence the probability 
cannot be expcx). 

Then, are the stable histories preferred? We have to compare the left-rolling (stable) histories 
and the right-rolling (unstable) histories. As we discussed above, our numerical results and the 
understanding of the behavior of the histories contributing to the saddle point approximation that 
we have developed in the previous sections show that the probabilities for a classical universe with 
stabilized, and one with non-stabilized couplings should be of similar magnitude. At very least, it 
appears that the stabilized universe would not be exponentially suppresed, unless the potentially 
would somehow be extremely ("exponentially") asymetric. 

In conclusion, we can say on the dilaton stabilization problem using the no-boundary measure 
that: 

1. We disagree with the assertion that the probability of a universe corresponding to the runaway 
solution is 1. 

2. Rather, under the approximation we are working in, probabilities for stabilization and de- 
stabilization are of similar order. 

Therefore, the no-boundary proposal seems to be able to contribute to the solution of the dilaton 
stabilization problem. 

5 Discussion 

In this paper, we investigated the no-boundary measure in scalar-tensor gravity in the context of 
Euclidean quantum cosmology. Especially, we were interested in trying to explain why we see a 
non- vanishing gravitational coupling. This is related to the dilaton stabilization problem, to explain 
why the dilaton field is located in a stable vacuum of a potential. 

To get a handle on this problem, we worked with a Brans-Dicke field with various potentials. 
As was found in previous work, we have seen that it is important to not only consider real field 
configurations in the saddle point approximation to the wave function, but also complex ones, the 
fuzzy instantons. In our case, we see that fuzzy instantons contributing to the probability for 
classical universes are allowed in the following two regions in a potential: With the scalar at a local 
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minimum or near a local m.axim.um. The instantons at the m.inimum. is in fact not fuzzy, but they 
are of measure zero as compared to the other fuzzy instantons. Therefore, even though the former 
have large negative Euclidean action, their probability will be zero. 

Thus, histories contributing to the probability for a classical universe should slow down and 
spend time around a local maximum, experience false- vacuum inflation, and then turn to Lorentzian 
time. This is indeed a new classicalization mechanism that is different from the previous work of 
Hartle, Hawking, and Hertog. 

The next question is whether we will roll-down to left or to right. In at the level of precision 
and approximation that we are working at, we cannot decide. Probabilities of two possibilities 
have similar order, unless the potential is extraordinarily asymmetric. Therefore, the no-boundary 
measure can partly explain the stabilization of some coupling constants of nature. But it docs not 
seem to assign probability zero universes in which they run away. This is perhaps a point where 
anthropic reasoning may be employed. 

There are some issues that we have not studied, but which may nevertheless be relevant to the 
questions we considered: 

• In this work, we did not include volume weighting, since we did not include an inflaton field 
(Of course, the dilaton field can role as an infiaton field, but it may not necessarily be true). 
If we include the inflaton field, the classicalization process could be changed, since there is 
another field that induces inflation. The inclusion of volume weighting and/or an additional 
scalar may change the results. 

• If there is a correction term which affects histories near a local minimum so that it breaks 
symmetry around the local minimum, then the region near the local minimum could admit 
a continuous spectrum of histories. Then, if vacuum energy of such a local minimum is 
sufficiently smaller than other positions, it may fully explain the stabilization of the dilaton 
field. 

• From our numerical results it seems that our intuition, coming from the dynamics of real fields, 
serves us well for the qualitative understanding of the dynamics of the fuzzy instantons. It 
should however be kept in mind that these are complex field configurations. In particular, 
they are sensitive to the analytic continuation of the potential. It is thus not inconceivable 
that different potentials, with a similar shape on the real sub-sector of the theory, will lead 
to different results. 

In addition, although our study is on scalar-tensor gravity, our conclusion is qualitatively relevant 
also for the no-boundary measure of Einstein gravity minimally coupled to a scalar field. Due to 
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their choice of the potential, Hartle, Hawking and Hertog [51 [7] only considered slow-roll inflation, 
concluded that the bottom-up probabilities do not favor larger amounts of inflation, and then 
appealed to volume weighting to deal with this problem |3j. Our work shows, however, that there 
is a second mechanism to generate classical histories, provided the potential has a more complex 
shape. They can be formed around the local maximum via false vacuum inflation. In this case, the 
bottom-up probability does not disfavor the top of the hill of the potential, and it may be possible 
to explain large amounts of inflation from bottom up probabilities in this way. 
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